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To study the optimum design factors and the optimal thermo-economical performances an optimal 
performance analysis of a solar driven heat engine system at maximum thermo-economical objective 
function conditions is done numerically. In the present investigation, thermodynamic analysis and 
NSGAII algorithm are employed to optimize objective function associated to the power output, thermal 
efficiency for a Solar driven engine system. Three decision-making procedures are applied to optimized 
answers from the results. The error through investigation is shown using error analysis. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Performance optimization studies of heat engines using finite¬ 
time thermodynamics were begun by Chambadal [1] and Novikov 
[2] and were continued by Curzon and Ahlborn [3]. At first step, the 
performance of an endoreversible Carnot heat engine at maximum 
power output are scrutinized by Curzon and Ahlborn [3 . Over the 
recent decade, great number of optimization investigations for 
heat engines based on endoreversible and irreversible models have 
been done owing to finite-time and finite-size constraints under 
various heat transfer conditions, commonly linear and non-linear 
ones. The eager authors could refer to a literature survey compiled 
by Bejan [4, Chen et al. [5] and Durmayaz et al. [6 . The power out¬ 
put, thermal efficiency, entropy generation, the ecological benefit 
and thermo-economic objectives are often selected as optimization 
factors. Great attention has been drawn to optimization of the so¬ 
lar-driven heat engines based on various performance criteria [7- 
21 where efficiency, power output and power density are often se¬ 
lected as the objective functions for optimization. The optimum 
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operation conditions of an endoreversible and irreversible solar 
driven heat engine was studied by Sahin [17,18] and Sogut and 
Durmayaz, respectively [20]. 

Sahin et al. [22], employed this method to study the thermoeco¬ 
nomics of an endoreversible heat engine in terms of the maximiza¬ 
tion of a profit function defined as the quotient of the power output 
and the annual investment cost. Recently, Barranco-Jimenez et al. 
[23], analyzed the optimum operation conditions of an endorevers¬ 
ible heat engine with different heat transfer laws at the thermal 
couplings but operating under maximum ecological function con¬ 
ditions. Furthermore, the thermoeconomic optimum operation 
conditions of a solar-driven heat engine have been investigated 
by Barranco-Jimenez et al. [24]. 

In these studies, they considered three regimes of performance: 
(i) The maximum power regime (MPR) [25-27], (ii) the maximum 
efficient power [28,29 , and (iii) the maximum ecological function 
regime (MER) [30,31]. 

In this work, we study the thermoeconomics of an irreversible 
heat engine by considering losses due to heat transfer across finite 
time temperature differences [32-35 , heat leakage between ther¬ 
mal reservoirs [36-44] and internal irreversibilities [45-47] in 
terms of a parameter which comes from the Clausius inequality. 
In this work, two regimes of performance are considered: The 
maximum power regime, and the ecological function regime. In 
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Nomenclature 



A 

area 

e 

T h 

T h 

a, b 

investment cost proportionality coefficient 



c, 

annual investment cost 

Subscrints 

Fij 

objectives for solutions of the Pareto frontier 

H 

related to high temperature thermal reservoir 

F n . 

y 

non-dimensionalized objective 

L 

related to low temperature thermal reservoir 

Pm 

objective function associated to the power output 

h 

hot side 

i 

relative investment cost 

c 

cold side 

W 

output power 

R 

non-endoreversibility parameter 

Q 

heat transfer rate 

LK 

leakage 

T 

temperature 

HT 

total heat rate 

T 

1 a 

the temperature of the fluid surrounding the body 



u 

heat transfer coefficient 

Acronyms 



NSGA 

Nondominated Sorting Genetic Algorithm 

Greek letter 

MPR 

maximum power regime 

a 

thermal conductance 

MER 

maximum ecological regime 

p 

ratio of the heat conductance parameter 

EA 

evolutionary algorithm 

T 

heat sources temperature ratio 

MOEAs 

multi-objective evolutionary algorithms 

n 

efficiency 



y 

internal conductance 



( 

ratio of internal conductance 




addition, in the present work, thermo-economical analysis, con¬ 
ductive-convective and radiative terms are considered through a 
heat transfer law of the Dulong-Petit type [48,49 . 

Multi-objective optimization has supreme presence in different 
engineering problems such as skyline computation, vehicle routing 
issues and so forth [53-55 . Solution of the multi-objective optimi¬ 
zation issues is an acutely hard purpose which needs the contem¬ 
poraneous satisfaction of a number of different and sometimes 
conflicting objectives. Evolutionary algorithms (EA) were ex¬ 
panded and applied during the 18th century in an effort to stochas¬ 
tically solve problems of this generic class [56 . A reasonable 
solution to a multi-objective quandary is to inquiry a set of solu¬ 
tions, each of them satisfies the objectives at aplausible level with¬ 
out being overcome by any other solution [57]. Multi-objective 
optimization quandaries universally represent a possibly uncount¬ 
able collection of solutions so-called Pareto frontier, whose evalu¬ 
ated vectors illustrate the best feasible trade-offs in the objective 
function space. In this case, multi-objective optimization of differ¬ 
ent themodynamic and energy systems have been drawn attention 
by investigators these days [58-65]. In the present investigation, 
multi-objective optimization is conducted to achieve the optimal 
performance analysis of a solar driven heat engine system at max¬ 
imum ecological objective function conditions. In our study, we 
consider different heat transfer laws in the thermal couplings 
and several decision variables such as temperature ratio, heat 
source temperature, hot working fluid temperature, ratio of heat 
transfer areas (heat sink to heat source) are considered. The pre¬ 
sented analysis is great at designing and suitable design parame¬ 
ters can be obtained based on the analysis. The paper is 
organized as follows: In Section 2, we present the thermodynamic 
description of the system, we apply the first and second laws of 
thermodynamics to construct the objective functions, in Section 3 
we describe the multi-objective optimization using evolutionary 
algorithms, in Section 4, we describe our decision-making methods 
used in our work, in Section 5 we present our main results and fi¬ 
nally in Section 6, we present our conclusions. 

2. Theoretical model 

The considered irreversible solar-driven heat engine operates 
between a heat source of temperature T H (the collector) and a 


heat sink of temperature T L (cooling water), see Fig. la. Working 
fluid temperature which exchange heat with reservoir at temper¬ 
ature T h is cosidered as T h . Similarly, working fluid temperature 
which exchanges heat with reservoir at temperature T L is defined 
as T c . A T-S diagram of the model including heat leakage, finite 
time heat transfer and internal irreversibilities is also shown in 
Fig. lb. 

An attempt to describe combined conductive-convective and 
radiative cooling by a power-law relationship is given by the so- 
called Dulong-Petit law of cooling [23 , which is 

Q = UA(T a -T) n (1) 

where Q is the rate of heat loss per unit area from a body at temper¬ 
ature T, T a is the temperature of the fluid surrounding the body, and 
n is an exponent with a value between 1.1 and 1.6 [23 , Some 
authors have established that n = 5/4 based on studies made by Lor- 
entz and Langmuir. In the present paper, we use the DP law of cool¬ 
ing with n = 5/4. 

The rate of heat leakage Q LI < from the hot reservoir at tempera¬ 
ture T h to the coldreservoir at temperature T L with thermal con¬ 
ductance y is given by [48] 

Qlk = y(T H - r L f = iU H A H (T H - T L f (2) 

where y is the internal conductance of the heat engine and £ de¬ 
notes the ratio of the internal conductance with respect to the 
hot-side convection heat transfer coefficient and heat transfer area, 
that is, £ = [47 . The 5/4 exponent is usual in a Dulong-Petit 

heat transfer law [48]. The rate of heat flow from the hot source 
to the heat engine is given by, 

0.H = U„A„(T H - T h )i (3) 

where U H is transfer coefficient and A H is the heat transfer area of 
the hot side heatexchanger. Eqs. (2) and (3) are of the Dulong-Petit 
type, including conductive-convective and radiative effects [48] 
rather than only use a Stefan-Boltzmann law. Note that, Stefan- 
Boltzmann law is proper for spaceapplications [31 , where the effect 
of atmospheric gases is not present. On the other hand, a Newtonian 
heat transfer is presumed as the main mode of heat transfer to the 
low temperature reservoir, so the heat flow rate from the heat en¬ 
gine to the cold reservoir, Ql, can be calculated as, 
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Fig. 1. Schematic diagram of the irreversible heat engine and its T-S diagram. 


Ql = U l Al(T c - T L f (4) 

where U L is the cold side heat transfer coefficient and A L is the heat 
transfer area of the cold side heatexchanger. Then, the total heat 
rate Qht transferred from the hot reservoir is, 

Qht = Q.h + Qlk = U h A h (T h — Th) 4 + £UhA h (T h — Ti) 4 (5) 

and the total heat rate Qlt transferred to the cold reservoir is, 

Qlt = Ql + Qlk = U l A l (T c - T l ) + f U h A h (T h - T L f 4 (6) 

Applying the first law of thermodynamics, the power output is 
given by, 

W = Qht - Qlt = U H A H (T H - T h f - U L A L (T C - T L ) (7) 

Applying the second law of thermodynamic to the irreversible 
part of the model we get, 


dQ Qh Ql 


<0 


T T h T c 
One can rewrite the inequality in Eq. (7) as, 



Qh _ R Qi 


( 9 ) 


where R is the so-called non-endoreversibility parameter [50,51 j. 
Substituting Eqs. (3) and (4) into Eq. (9), a relationship between T c 
and T h is obtained as, 


_ _ RAr _ 

Tl (ra r - 



where 6 = ?£, A R and p are the ratios of the heat transfer areas and 
the heat conductance parameter respectively, and are defined as, 



Al 

Ah 



U h i 

P = JJ^ T h ( 12 ) 

These two parameters can be taken as design parameters. The 
thermal efficiency of the irreversible heat engine is, 


_-j_ Qlt _ Qh — Ql 

Qht Qh + Qlk 



In thermoeconomic analysis of power plant models, an objec¬ 
tive function can be specified in terms of a characteristic function 
(power output [18,26,52], ecological function [24,30,49], etc.) and 
the cost involved in the performance of the system. 

In this context, De Vos [26] studied the thermoeconomics of a 
Novikov power plant model in terms of the maximization of an 
objective function expressed as the quotient of the power output 
and the operating costs of the plant. 

In that paper, De Vos considered a function of costs with two 
parts: (i) The cost of the investment which is assumed as propor¬ 
tional to the size of the plant and (ii) the cost of the fuel consump¬ 
tion which is assumed to be proportional to the heat input in the 
Novikov model, similarly, Sahin and Kodal base on Curzon and Ahl- 
born [25] model, formulated a thermoeconomic analysis in terms 
of an objective function. Their function defined as power output 
per unit total cost including both the investment and fuel costs 
[18], but assuming that the size of the plant can be taken as pro¬ 
portional to the total heat transfer area, instead of the maximum 
heat input which previously considered by De Vos [26]. 

Following the Sahin et al. procedure [52], the objective function 
has been defined as the power output per unit investment cost due 
to a solar driven heat engine does not consume fossil fuels. To opti¬ 
mize power output per unit total cost, the objective function is de¬ 
fined by [52] as follows: 



where C, represents annual investment cost of plant. C, is assumed 
to be proportional to the size of the plant. Also, the size of the plant 
can be proportional to the total heat transfer area. Consequently, 
the annual investment cost of the system can be written as [52], 

C; = clAh + bAi (15) 


where the investment cost proportionality coefficients for the hot 
and cold sides a and b respectively are equal to the capital recovery 
factor times investment cost per unit heat transfer area, and their 
dimensions are ncuKyear m 2 ), ncu being the national current unity. 
By using Eqs. (3), (7), (10), (14), and (15), we get a normalized 
expression for the objective function associated to the power output 
given by [52], 

w (j L -h 


r M = r = 



(16) 
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where t = ^ and the parameter/, is the relative investment cost of 
the hot size heat exchanger defined as [52 , 


/ 


a 

a + b 



And, by using Eqs. (2)—(4), the thermal efficiency, rj } of the irre¬ 
versible heat engine can be expressed by, 


xP(\-ey-A R fe-\) 

11 =- ? (18) 

t/?(1 -0)’ + £t 5 


3. Multi-objective optimization with evolutionary algorithms 

3.1. General concepts of the multi-objective optimization 

Imagine a decision-maker that would like to optimize objec¬ 
tives which are non-commensurable. The decision-maker does 
not know the priority of the objectives. In addition, all objectives 
are of the minimization kind. For a minimization type of objectives, 
it can be altered to the maximization type by multiplying negative 
one. A minimization (optimization) multi-objective decision prob¬ 
lem with K objectives is defined as follows. 

Givenan non-dimensional decision variable vector x = {xi,- 
..., x n ] in the solution area X, find a vector x* which minimizes a 
collection of I< objective functions/(x) = {/j(x*),... ,/<(x*)}. The solu¬ 
tion space X is restricted by a number of constraints, such as gf- 
x*) = bj for j = 1,..., m, and bounds on the decision variables [57]. 
Totally, it is not possible to minimize all the I< objective functions 
simultaneously with one solution vector X which means in real- 
world problems, objectives underconsideration conflict each other. 
Thus, optimizing X with concerning a uni-objective mostly leads to 
unreasonable consequence with respect to the other objectives. 
Therefore a novel implication, so-called the “Pareto optimum solu¬ 
tion”, is applied in multi-objective optimization quandaries. A 
practical solution X is called “Pareto optimal” if there is no other 
plausible solution Y that prevails solution X. Therefore, a feasible 
solution Y is entitled to overcome another conceivable solution X, 
if and only if, /(Y) ^ /(X) for i = 1,..., k and f(Y) </(X) for at least 
one objective function which explains that a possible vector X is 
denominated Pareto optimal if there is no other possible solution 
Y that would decrease some objective functions without conclud¬ 
ing a simultaneous increase in minimum one other objective func¬ 
tion. The collection of all plausible non-dominated solutions inXis 
leaded to the “Pareto optimal set”, and for a certain Pareto optimal 



collection, the corresponding objective function values in the 
objective space are denominated the “Pareto optimal frontier” 

[57]. 

Fig. 2 represents the objective functions space for an optimiza¬ 
tion problem with two objective functions f\ and / 2 . The feasible 
are is insight and infeasible area is outsight of the curve. Each point 
in the feasible area denotes a solution. It is obvious that the values 
of functions/i and f 2 for point M are lower than the corresponding 
values of point j. Therefore point M overcomes point J or point M is 
better than point J. In the same way points L and N dominate M. But 
points I and I< neither overcome M nor are overcome by M. Thus, 
only the points that are placed in the left-down parts of M would 
overcome M. There is not any point in the feasible space placed 
in the left-down part of point R. Thus R is not overcome by the 
other points in the feasible space and hence R is a Pareto optimal. 
This is true for points P, A, R, E, T, 0 and all of the other points lo¬ 
cated at the bold curve interpreted as “The Pareto Optimal Fron¬ 
tier” in Fig. 2. In the feasible space, the minimum value of /i is 
for P and the minimum value of f 2 is for point 0. Thus P and 0 
are the solutions for one objective optimization issues whose 
objective functions are /i and / 2 , respectively. The other points on 
the Pareto frontier are also optimal, but a decision-making proce¬ 
dure is needed to find which of them should be choosen. 

3.2. Optimization via EA 

In the present investigation, the Pareto frontier is obtained 
using the Genetic Algorithm (GA) which is identified as a field of 
evolutionary algorithm, genetic algorithms were evolved by John 
Holland in the 1960s to import the mechanisms of natural adapta¬ 
tion into computer algorithms and numerical optimization [59]. 
They are employed as a computer simulation in a way that a pop¬ 
ulation of abstract representations (called chromosomes or the 
genotype of the genome) of candidate solutions (called individuals, 
creatures, or pheno-types) to an optimization problem evolves to¬ 
ward better solutions. It commences from a population of ran¬ 
domly produced exclusives and occurs in generations. In each 
generation, the fitness of every individual is measured; multiple 
individuals are stochastically opted from the current population, 
and altered to form a new population. The new group is employed 
in the next iteration. Generally, the algorithm finishes when either 
a maximum number of generations have been produced, or a sat¬ 
isfactory fitness level has been accomplished for the population. 
If the algorithm has ended owing to a maximum number of gener¬ 
ations, a satisfactory solution might have been reached or might 
not have been gained. In genetic algorithms, a candidate solution 
to a problem is named a chromosome, and the evolutionary viabil¬ 
ity of each chromosome is denoted by a fitness function. This 
method can be applied to optimize non-linear problems [57,60 . 

Multi-objective evolutionary algorithms (MOEAs) have been 
evolved in the recent decades by a lot of tests on complex mathe¬ 
matical quandaries. On real-world it has shown that they can omit 



Fig. 2. Schematic of the objectives space. 


Fig. 3. Scheme for the multi-objective evolutionary algorithm used in the present 
study. 
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the obstacles of classical methods [57,60]. The structure of the 
MOEA exploited in this investigation is depicted in Fig. 3 [59]. 
The real values of decision variables are applied instead of their 
binary coded. 

3.3. Objective functions, decision variables and constraints 

Two objective functions for this study are the thermal efficieny 
( rj ) and the objective function associated to the power output (F M ) 
respectively. 

In this paper five decision variables have been considered as 
follows: 

6: temprature ratio (e = ^, 

A r : ratio of the heat transfer areas, 
t: temperature ratio (j£), 

/: the relative investment cost of the hot size heat exchanger, 
p: ratio of the heat conductance parameter. 

The objective functions with respect to following constraints 


have been solved: 


0.4 < 9 < 0.9 

(19) 

2 ^ A r ^ 4 

(20) 

3 ^ T ^ 5 

(21) 

0.5 0.9 

(22) 

0.5 ^ p^ 1.5 

(23) 


4. Decision-making in the multi-objective optimization 

In multi-objective optimization a process of decision-making 
for selection of the final optimal solution from available solutions 
is needed. There are a number of methods for decision-making pro¬ 
cess in decision problems. These methods can be utilized for selec¬ 
tion of a final optimum solution from the Pareto frontier. Since, 
dimension of various objectives in a multi-objective optimization 
problem might be different, therefore, before any decision, dimen¬ 
sions and scales of objective space should be unified. In this regard, 
objectives vectors should be non-dimensioned before the decision¬ 
making process. There are various methods of non-dimensionaliza- 
tion employed in decision making including Euclidean non-dimen- 
sionalization, and fuzzy non-dimensionalization. 

4.1. Non-dimensionalization methods 


4.1.1. Euclidean non-dimensionalization 

The matrix of objectives for various solutions of the Pareto fron¬ 
tier is denoted by Fy where i is an index for each solution on the 
Pareto frontier and j is an index for each objective in objective 
space. In this method, a non-dimensionalized objective, FjJ, is de¬ 
fined as, 



J = For minimizing and maximizing objectives 



4.1.2. Fuzzy non-dimensionalization 

In this method, a non-dimensioned objective, Fj|, is defined as, 



Fy - min(Fg) 

max(Fy) - min(Fy) 


(25a) 



max(Fy) - Fjj 
max(Fy) - min(Fy) 


For minimizing objectives 


(25b) 


In this paper, most familiar and usual type of decision-making 
processes including the fuzzy Bellman-Zadeh, LINMAP and TOPSIS 
method are applied in parallel and the final optimal solution is se¬ 
lected using these three methods. The Bellman-Zadeh method uti¬ 
lizes the fuzzy non-dimensionalization while the LINMAP and 
TOPSIS methods utilize Euclidean non-dimensionalization. The fol¬ 
lowing sections are described these decision-making algorithms. 


4.2. Decision making methods 

4.2.1. Bellman-Zadeh decision-making method 

A final decision is defined by the Bellman and Zadeh model as 
the intersection of all fuzzy criteria and constraints and is repre¬ 
sented by its membership function. In this method, a matrix of 
membership function is defined. Column of this matrix includes 
the fuzzy membership function for each objective. Rows of the 
membership matrix reflect magnitudes of the membership func¬ 
tions for each solution of the Pareto frontier. Therefore, the number 
of columns in the membership function matrix is equal to the 
number of objectives and the number of rows is equal to the num¬ 
ber of solutions located on the Pareto frontier. Detail for method of 
definition for the membership function has been presented in Ref. 
[66]. The fuzzy membership function for each solution is set at the 
minimum value of the membership functions of all objectives for 
the proposed solution. Therefore, a vector of membership is 
obtained in which its elements are the minimum membership 
function of objectives at each solution. Finally, the maximum of 
these minimum, i.e. an element of the membership vector with 




For maximizing objectives 


Fig. 4. (a and b) Pareto optimal frontier in objectives’ space. 
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Population Population 

Fig. 5. Scatter distribution of 9 with population in Pareto frontier. Fig. 8. Scatter distribution of p with population in Pareto frontier. 



Population 



Fig. 6. Scatter distribution of t with population in Pareto frontier. 


Fig. 9. Scatter distribution of / with population in Pareto frontier. 



0 50 100 150 200 250 300 350 400 450 500 


Population 

Fig. 7. Scatter distribution of A R with population in Pareto frontier. 


maximum values of membership function is selected as a final 
solution. In other word, in the Bellman-Zadeh fuzzy decision mak¬ 
ing a maximum of minimums (minimum membership function of 
objectives) is selected as the final optimal solution. In Ref. [66] 
more detail is presented. 


impossible to have each objective in its optimal condition that it 
can obtain in a single-objective optimization. Consequently, the 
ideal point is not located on the Pareto frontier. In the LINMAP 
method, after Euclidean non-dimensionalization of all objectives, 
the spatial deviation of each solution on the Pareto frontier from 
the ideal point is determined. Finally, a solution with a minimum 
spatial distance from the ideal point is selected as a desired final 
optimal solution. More details of the LINMAP method have been 
presented in Ref. [66 . 


42.3. TOPSIS decision-making method 

In this method, beside the ideal point, a “non-ideal point” is de¬ 
fined. The non-ideal point is the ordinate in objective space in 
which each objective has its worst value. Therefore, beside the 
solution spatial distance from the ideal point, the solution spatial 
distance from the non-ideal is utilized as a criterion for selection 
of the final solution. In the TOPSIS method a solution with a max¬ 
imum distance from the non-ideal point and minimum distance 
from the ideal point is selected as a desired final optimal solution. 
More details of the TOPSIS method have been given in Ref. [66]. 


5. Results and discussion 


4.2.2. LINMAP decision-making method 

An ideal point on the Pareto frontier is a point in which each 
objective is optimized regardless to the satisfaction of other 
objectives. It is clear that in the multi-objective optimization it is 


The objective function associated to the power outputand ther¬ 
mal efficiency of the solar driven engine are maximized simulta¬ 
neously using the multi-objective optimization method which 
works based on the NSGA-II algorithm. 
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Table 1 

Decision making of multi-objective optimal solutions. 


Decision maker 

Decision variables 





Objective functions 


e 

Ar 

T 

/ 

P 

n 

Fm 

TOPSIS 

0.671637 

3.997406 

4.99973 

0.89996 

1.326286 

0.5356731 

0.52153751 

LINMAP 

0.679862 

3.990805 

4.999976 

0.899967 

1.336159 

0.542754 

0.513691706 

Bellman-Zadeh 

0.722324 

3.992746 

4.99979 

0.899942 

1.472895 

0.5735338 

0.46678106 


In this regard, optimization is performed with objective func¬ 
tions that are expressed by Eqs. (16) and (18) constraints which 
are expressed with Eqs. (19)-(23). 

Design parameters (decision variables) of optimization are the 
temperature ratio (o = ^, ratio of the heat transfer areas, temper¬ 
ature ratio (r = 17 ), the relative investment cost of the hot size 
heat exchanger, ratio of the heat conductance parameter. 

Pareto optimal frontier for two objective functions, objective 
function associated to the power output, thermal efficiency of Solar 
driven engines are represented in Fig. 4a. The chosen points based 
on decision making methods are shown, too. The fitted curve equa¬ 
tion of the Pareto optimal frontier is obtained and illustrated in 
Fig. 4b, which can be applied to gain system thermal efficiency, 
as following: 

„ = (~2163F M + 1389) (25c) 

(F 2 m - 3367 F„ + 2246] 

R 2 = 0.9982 
0.46 ^ Fm ^ 0.62 
0.25 ^ rj ^ 0.55 

Figs. 5-9 demonstrate the scattered distribution of the design 
variables. These plots can provide better vision of variation of the 
design variables from the Pareto frontier. According to this distri¬ 
bution, it is found that heat sources temperature ratio and relative 
investment cost reach their maximum values. It displays that in¬ 
crease in these two design parameters result in improvement in 
both objective functions. Consequently, their maximum values 
are chosen. 

The optimal results for decision parameters and objective func¬ 
tions using LINMAP, TOPSIS, Bellman-Zadeh decision-making 
methods are summarized in Table 1. 

Maximum and mean error analyses for objective functions are 
categorized in Table 2. The analysis is done for all mentioned deci¬ 
sion-making methods. 

First row of Table 2 express maximum absolute percentage er¬ 
ror of three decision making methods and second row of this table 
report mean absolute percentage error of each decision making 
methods. In both cases errors are less than 1%. 

The mean absolute percentage error was calculated in the error 
analysis of implementing methods. To assess this goal, 30 runs of 
computer code were performed and a final solution selected by 
fuzzy Bellman-Zadeh, LINMAP and TOPSIS decision making meth¬ 
ods were obtained for each run. Absolute percentage error of solu¬ 
tions obtained in each run was calculated comparing to the best 
results introduced as the final solution in this paper. 


Table 2 

Error analysis based on the mean absolute percent error (MAPE) method. 


Errors of decision makers 

LINMAP 

TOPSIS 

Bellman-Zadeh 


n 

Fm 

n f m 

n F m 

Max error % 

0.05 

0.29 

0.39 0.23 

0.19 0.10 

Average error % 

0.03 

0.20 

0.15 0.11 

0.11 0.06 


6. Conclusions 

Throughout current research work enormous efforts were made 
to figure out the optimum condition of the two addressed objective 
functions such as thermal efficiency and the objective function 
associated to the power output for a solar driven heat engine. In 
developed multi-objective optimization approach five different 
variables including temperature ratio (0 = ^, temperature ratio 
(t = ^, ratio of the heat transfer areas (A R ), ratio of the heat con¬ 
ductance parameter (/?), the relative investment cost of the hot size 
heat exchanger (J) were considered as decision variables. Due to 
obtained solutions from different decision making approaches, TO- 
SIS decision making has higher maximum error percent than other 
two addressed decision making approaches such as LINMAP and 
Fuzzy approaches. The results of the presented thermodynamic 
model for solar driven engine implementing NSGAII algorithm 
are acceptable. 
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